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CHARTS AND APPROXIMATE FORMULAS FOR THE ESTIMATION ОЕ AEROELASTIC EFFECTS 
ON THE LOADING OF SWEPT AND UNSWEPT WINGS ! 


“Ву Franxuin W. DrEDERICH and KxwNNETH А. Foss 


SUMMARY 


Charls and approximate formulas are presented for the 
estimation of aeroelastic effects on Ше spanwise lift distri- 
bution, lift-curve slope, aerodynamic center, and damping 
in roll of swept and unswept wings at subsonic and supersonic 
speeds. Some design considerations brought out by the results 
of this report are discussed. 


INTRODUCTION 


A knowledge of the spanwise lift distribution and of some 


of the aerodynamic parameters associated with it is required - 


for the design of a wing structure. Under certain conditions, 
such as high dynamic pressures, thin wings, swept wings, 
or wings designed for low wing loadings, the spanwise lift 
distribution may be affected to a significant extent by 
aeroelastic effects, because a wing which carries a certain lift 
necessarily deforms under that lift. If the angles of attack 
along the span are changed as a result of this deformation, the 
lift carried by the wing is changed as well; in turn, this 
change in lift causes a change in the deformation of the wing 
and hence another change in lift, and so on, until an equilib- 
rium condition is reached. The changes in the magnitude 
and the distribution of the lift are reflected in changes of the 
wing lift-curve slope, the wing bending and rolling moments, 
the spanwise center of pressure of the lift, and, on a swept 
wing, the longitudinal center of pressure. 

Inasmuch as the lift produced by a given change in angle 
of attack is proportional to the dynamic pressure, the various 
acroelastic effects tend to increase with dynamic pressure. 
In fact, for certain wings а sufficiently large dynamic pressure 
may produce a condition of instability in which the change in 
lift caused by deformation is greater than the amount of 
lift required to produce the deformation, so that a given 
deformation will tend to increase until the structure fais. 
This phenomenon is aeroelastic divergence; since it involves 
only torsional deformations in the case of unswept wings, itis 
often referred to as torsional divergence. 

Several methods are available for calculating these effects 
(ref. 1, for instance), but since these effects depend on the 
structural characteristics of the wing, which are not ас- 


curately known in advance of its design, the relatively large 
amount of time required for even the most efficient of these 
methods militates against their use in connection with 
preliminary design calculations. A need exists, therefore, 
for means of estimating some of the more important aero- 
elastic effects on the spanwise lift distribution quickly and 
with an accuracy that is sufficient for preliminary design 
purposes. 

Charts and approximate formulas are presented in this 
report for estimating the changes in spanwise lift distribution, 
lift-curve slope, wing rollmg-moment coefficient, spanwise 
center of pressure, and aerodynamic center occasioned by 
&eroelastic action of swept and unswept wings at subsonic 
and supersonic speeds. Also included are summary charts 
which indicate whether the various aeroelastic phenomena 
considered are likely to affect any given design. By means 
of these charts the conventional procedure of designing a 
wing on the basis of certain strength criteria, checking it for 
aeroelastic phenomena, and then reinforcing it, when neces- 
sary, to meet the stiffness requirements imposed by these 
phenomena can often be simplified greatly, masmuch as 
the effect of some of these phenomena can be estimated in 
advance of design. 

The use of the charts is described in the section headed 
“Calculation of the Various Aeroelastic Phenomena,” and 
some considerations involved in the selection of the aero- 
dynamic, structural, and geometric parameters are discussed 
in some detail in the section headed “Selection of Param- 
eters.” These two sections, аз well as the sections headed 
"*[lustrative Example” and “Preliminary Survey of Aero- 
elastic Behavior," are likely to prove of greatest interest 
аф а first reading of this report. The various parts of the: 
section headed “Discussion” are concerned with the limita- 
tions of the charts, with the light they shed on such practical 
design problems as the relative significance of strength and 
stiffness as design criteria, with efficient ways of stiffening 
a wing that is strong but not stiff enough, and with the 
achievement of aeroisoclinic conditions. 

A. brief description of the calculations (based on refs. 1 
and 2) used in preparing the charts is contained in the 
appendixes. 


1 Previously released as МА ОА ТМ 2608, “Charts and Approximato Formules for the Estimation of Aeroelastic Effects on Ше Loading of Swept and Unswept Wings," by Franklin W. 


Diederich and Kenneth А, Foss, 1952, 
321000-50----47 
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SYMBOLS 


aspect ratio, 03/9 

swept-span aspect ratio, АТсозА 

cross-sectional area of the (assumed) single 
torsion cell, sq in. 

distance from leading edge to section aerody- 
namic center, fraction of chord 

distance from leading edge of mean aerodynamic 
chord to wing aerodynamic center, fraction of 
mean aerodynamic chord . 

wing span, in. 

wing span less width of fuselage, em in. 

wing-root bending-moment coefficient, 4 М. Ј985 

lift coefficient of wings alone (exclusive of fuse- 
lage), Le/qS 

wing lift-curve slope per radian 

effective lift-curve slope per radian (defined in 
eqs. 2 and 4) 


го -moment coefficient on both wings alone Қ 


(exclusive of fuselage), Rolling moment/qSb 
wing-root twisting-moment coefficient, 27;,/qSc, 
chord (measured perpendicular to elastic axis), 

in. | 
ote С: . Ша 


2 ІП. 


average chord, 


section lift-curve да рег radian 

mean aerodynamic chord (parallel to plane of 
symmetry), in. 

Young's modulus of'elasticity, Ib/sq in. 

distance from leading edge to elastic axis, frac- 
tion of chord 

dimensionless moment arm of the section lift 
about the elastic axis, e—a 

effective or average dimensionless moment arm 

allowable bending stress, Ib/sq in. 

root-stiffness function 

allowable shear stress, Ib/sq in. 

structural weight function 

dimensionless parameters used in approximate 
formülas for angle of attack due to aero- 
elastic deformation 

dimensionless functions -of the distance along 
the span used in approximate formulas for 


angle of attack due to aeroelastic deformation : 


modulus of rigidity, lb/sq in. 

wing thickness, in. 

section bending moment of inertia, in.* 

section moment of inertia in torsion, in. 

dimensionless parameters used in approximate 
formulas for dimensionless dynamic pres- 
sures at divergence 


dimensionless sweep param eter, - e 
Eu с; (ЕТ), 


lift of both wings alone O pi fuselage), Ib 
lift per unit distance along span, Бо. 


& (87), eua 


в < < < 5-24. 


ще 


>. 
“ 


>> pb nR 


bending moment about an axis perpendicular 
to elastic axis, in-lb 
free-stream Mach number 
design load factor 
rolling angular velocity, radians/sec 
dynamic pressure, lb/sq ft 
dimensionless dynamic pressure, . 
Cr. €1Cr 81" соз А 


144 (ӨЛ, 


dimensionless dynamic pressure, 
4 Сі, lrs” sin А 
144. (ЕЛ), 
total wing агеа, sq in. 
distance along elastic axis measured from wing 
root, in. 
dimensionless distance along elastic axis, 8/8, 
distance from root to center of pressure of lift 
along elastic axis, in. 


` dimensionless distance from root to center of 


pressure of lift, 8/8, 
accumulated torque about elastic axis, in-]b 
thickness of most highly stressed element of 
gkin, in. 
thickness of equivalent skin which includes the 
material in stringers and spar flanges, in. 
distributed torque due to inertia loading, in-]b/in. 
free-stream velocity, ft/sec 
design gross weight of airplane, Ib 
weight of primary structure of both wings, Ъ 
weight of both wings кын of fuselage, lb 
width of fuselage, in. 


| weight of primary ОТЕ УЕ structure per 


unit distance along span, Ib/in. 
lateral coordinate, in. 


dimensionless lateral coordinate, +2- тр 


lateral distance to center of pressure 
angle of attack in a plane parallel to plane of 


symmetry, radians 


‚ angle of local dihedral, radians; or spanwise 


slope of normal displacement of elastic axis 
density of the material of the primary structure 
(or an equivalent density in the case of sand- 
wich construction), Ib/cu in. 
lateral distance measured from wing root, 


„а 
y ә? In. 
dimensionless lateral distance, <= 


factors defined in table 1 

factors defined in equations (15а) and (15b) 
ratio of lift-curve slopes, Cz, 1 Cr, 

angle of sweepback at elastic axis 

wing taper ratio, с/с, 

free-air density, slugs/cu ft 
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e - angle of structural twist in planes perpendicular 
to elastic axis, radians 

7) tip stiffness ratio, (ЕГ), (ЕГ), 

шурут. dimensionless parameters used in approximate 


formulas for lift, root bending moment, and 
root twisting moment 


Subscripts: 

с8 constant stress 

D at divergence 

€ effective 

8 geometric (due to airplane attitude or built-in 
twist) 

a inertia 

0 rigid wing (for 9=0) 

r at wing root 

8 structural (due to structural or aeroelastic 
deformation) 

і at wing tip 

Superscrip ts: 

M due to bending moment 

Т due to torque 

Г due to root bending 

e due to root twist 


USE OF THE CHARTS AND APPROXIMATE FORMULAS 


SUMMARY OF METHOD AND ВСОРЕ OF CALCULATIONS ON WHICH THE 
CHARTS AND APPROXIMATE FORMULAS ARE BASED 


Although а detailed understanding of the method and 
scope of the calculations on which the charts of this report are 
based is not essential to the use of the charts, & brief account 
of these matters 13 given, primarily to эла in the appreciation 
of the limitations of the charts. Тһе method 18 described 
more fully in appendix А. 

Most of the calculations on which the charts are based 
were made by the method of reference 1, which consists in 
solving the differential equations descriptive of an elastically 
deformed wing under aerodynamic loading by numerical 
methods employing matrix techniques. Treated by this 


method were wings with four taper ratios А (1.0, 0.5, 0.2, and . 


0), two types of stiffness distributions (one proportional to 
the fourth power of the chord and one dictated by constant- 


stress considerations), and four values of a sweep parameter 


k at several values of the dynamic-pressure ratio rP Cal- 
D 


culated for each case were the dynamic pressure at divergence 
and the changes due to aeroelastic action in spanwise lift 
distribution, total wing lift, root bending moment, rolling 
moment, and spanwise center of pressure of the lift. For 
the wings of constant chord and constant stiffness, calcula- 
tions were also performed for six values of Е by а method 
which is an extension of that of reference 2 and consists in 
solving the differential equations exactly for these relatively 
simple cases. 

Some approximations have been made in the calculations 
concerning the aerodynamic induction effects, the root rota- 
tions, and the stiffness distributions, primarily in order to 
hold the number of variables considered in the analysis to a 
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minimum and to make the results more generally applicable. 

Aerodynamic induction effects at subsonic speeds are 
taken into account by an overall reduction of the strip- 
theory loading and, in the matrix calculations, by rounding 
off the strip-theory loading at the tip (see refs. 1 and 2); for 
supersonic speeds, strip theory is used with a small reduction 
at the tip in the matrix calculations. This approximation 
has made it unnecessary to consider explicitly the effects of 
aspect ratio, sweep, and Mach number on the rigid-wing lift 
distribution; the effects of these parameters on the total lift 
and on the aeroelastic increment to the lift distribution have 
been taken into account. - 

The rigid-body rotations imparted to a swept wing by its 
triangular root portion vary among different designs in а 
largely unpredictable manner. They have therefore been 
taken into account only by the use of an effective root, the 
selection of which in any given case is discussed briefly in a 
subsequent section. 

The spanwise distributions of the bending and torsional 
stiffnesses depend on the detailed design of the wing and 
cannot be generalized easily. The stiffness distributions 
used in the calculations of aeroelastic effects were obtained 
from the constant-stress concept outlined in appendix B, 
which constitutes an effort to relate the stiffness of a wing to 
its strength on the basis of the following assumptions: 

(1) The level of combined bending and torsional stresses 
is constant along the span. 

(2) The structure is designed for combined bending and 
torsional stresses in such a manner that the sum of the ratio 
of the actual to the allowable bending stress and the ratio of 
the actual to the allowable torsion stress is equal to unity 
when the margin of safety is zero. 

(3) The structure is of the thin-skin, stringer-reinforced 


. shell type and its main features do not vary along the span; 


for instance, the number of spars and their chordwise loca- 
tions are constant along the span. 

(4) At the design condition the spanwise distribution of 
the applied loading is proportional to the chord. 

Also used in the calculations were stiffness distributions 
which vary as the fourth power of the chord, as do those of 
solid wings and wings with geometrically similar cross sec- 
tions; as pointed out in a subsequent section, the results of 
these calculations can be used to estimate aeroelastic phe- 
nomena of some wings which have large cutouts or which 
for some other reason do not have stiffness distributions 
represented. fairly closely by those of the constant-stress 
type. 

All calculations are based on the assumptions that twisting 
is resisted primarily by the torsion cells of the wing structure 
and that the wing deformations can be estimated by means 
of the elementary theories of bending and torsion about an 
elastic axis. 

SELECTION OF PARAMETERS 

Geometric parameters.—The geometric parameters used 
in the analysis are defined in figure 1. The location of the 
effective root indicated in figure 1 is discussed in the section 
concerned with the structural parameters. 
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FIGURE 1.—Definitions of geometric parameters. 
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Aerodynamic parameters.— The aerodynamic parameters 
which enter the analysis are the wing lift-curve slope and 
the location of the aerodynamic center. Two lift-curve 
slopes are used at subsonic speeds: The wing lift-curve 
slope Cz, 18 used only in conjunction with additional lift 


distributions; for all other lift distributions—that is, those . 


due to built-in twist, due to roll, or due to aeroelastic twist— 
an effective lift-curve slope Cra, is used. Approximate 


values of these parameters are given for subcritical speeds 
by the relations 


А cos А 


C; ==, с; (1) 
А+2 == совА 
ай 
Сі =C; _ AcosA — (2) 


ав 


а с 
4-4 са Cos А 
Әт 


where сі, is the lift-curve slope of the section perpendicular 
to the quarter-chord line at а Mach number equal to M, cos А. 
Àn approximate value is given by 


РР 2т 
га J41— M соз?А 


Equation (1) is given in reference 3 and shown to be appli- 
cable both to incompressible and to subcritical compressible 
flow. Equation (2) is given in references 1 and 2 but without 
the term C; /2- in the denominator. This term is intro- 
duced into equation (2) in order to extend its applicability 


(3) 
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to compressible flows іп the same manner ав that employed 
for the coefficient of damping in roll in reference 3. (The 
lifting-surface corrections given in ref. 3 for the wing lift- 
curve slope and for the coefficient of damping in roll are not 
included in eqs. (1) and (2) because they are important 
primarily for wings of very low aspect ratio, to which the 
method of the present report is not applicable.) 

At supersonic speeds (more specifically, for supersonic 


leading and trailing edges) both lift-curve slopes are арргох- 
imately equal to the effective section lift-curve slope; that is 
= С = баш, 
where 
4 cos À 
= (4) 


С), ==-к======== 
“e «JMF сов? А—1 


The ratio of the lift-curve slopes С, and Ол, ів defined by 


k= | (5) 


so that for supersonic speeds к is equal to 1. 


The local aerodynamic centers are assumed to be at a 
constant fraction of the chord from the leading edge, so that 


' their distances from the leading edge (as fractions of the 


local chords) are all equal to the distance of the wing аего- 


‚ dynamic center from the leading edge’ of the mean aerody- 
namic chord (as a fraction of the mean aerodynamic chord). 


The moment arm e; is then given by the relation 
е=е—@, (6) 


The lift-curve slopes and the locations of Ше aerodynamic. 
center vary with the free-stream Mach number; hence the 
appropriate values must be used at each flight condition for 
which aeroelastic calculations are made. For airplanes 
designed to operate at subsonic speeds, only the highest 
Mach number attainable at the highest dynamic pressure is 
likely to be critical from aseroelastic considerations. Гог 
airplanes designed to operate at supersonic speeds no such 
general statement can be made; however, аба given altitude 
either the region of Mach numbers near the transition from 
the subsonic to the transonic regime or the highest attainable 
Mach number is likely to be critical as far as the aeroelastic 
phenomena considered in this report are concerned. (See- 
fig. 4 of ref. 1 and fig. 5 of ref. 2, for instance.) 

The airspeeds at which the various aeroelastic phenomena 
are of interest enter the calculations in the form of the cor- 
responding dynamic pressures. These dynamic pressures, 
in turn, are expressed in dimensionless form as 


Cr, ес 28, cos А 
а 


б Сас PORA 
чт OD, (7) 
ог 
С, с,8, gin A 
4 a (8) 
144 (ED, 


The parameter q* is useful in the analysis of unswept wings, 
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for which torsional deformations are predominant; the 
parameter 9 is primarily useful for highly swept wings, for 
which bending deformations are predominant. Іп general, 


the parameter q* is used in this report unless e; is zero. The 
ratio of these parameters 

i 4. 8; QD, 

k= 1* ao; (ED, tan A (9) 


is independent of the dynamic pressure and depends only on 
geometric and structural parameters. This ratio is very 
useful for analyzing the aeroelastic behavior of swept wings. 

Structural parameters.—For the purposes of an aeroelastic 
analysis, the wing structure is characterized by the location 
of the elastic axis and the magnitude and distribution of the 
bending and torsional stiffnesses (EJ and GJ), as well as by 
the magnitude of the rigid-body rotations imparted to the 
wing by its root. 

The elastic axis 18 usually defined аз the locus of points at 
which normal loads can be applied without causing the wings 
to twist, Such a locus does not geaerally exist for practical 
wings; however, for unswept wings without cutouts an axis 
can be determined which approximately satisfies this con- 
dition. Similarly, for swept wings without cutouts an elastic 
axis can be defined for the outboard part of the wing if the 
wing is considered to be clamped along some such line as the 
effective root shown in figure 1. In most aeroelastic calcula- 
tions the locus of shear centers for both swept and unswept 
Wings 18 assumed to be the elastic axis. If the structure has 
large cutouts which result in sudden changes in the stiff- 
nesses and in the shear center along the span, the charts of 
this report cannot be used except in a qualitative sense. 

The magnitude and the spanwise distributions of the 
bending and torsional stiffnesses enter aeroelastic calcula- 
tions by means of the charts and approximate formulas in 
different ways. The magnitudes, as characterized by the 
values of the stiffnesses at the effective root, have to be 
known in order to perform any calculations; the distributions 
are implicit in the charts. The root stiffnesses, if not known 
otherwise, can be estimated either from experience with 
similar designs, from the results of the constant-stress concept 
outlined in &ppendix B, or from & combination of the two. 

The required bending stiffness at the root (ЕГ), is propor- 
tional to the design load factor, the wing loading, the wing 
thickness ratio, the fourth power of the root chord, the 
square of the swept-span aspect ratio (А/совА), and the 
ratio of the modulus of elasticity to the allowable bending 
stress and depends on the taper ratio and on the detailed 
design of the wing (see appendix B). By means of this 
relation the bending stiffness of one wing сап be estimated 
from that of a reasonably similar wing. Or, with the con- 
stants of proportionality пъ and Ко given in table 1 and 
appendix B, respectively, which take into account some of 
the detailed design parameters as well as the taper ratio, 
the stiffness can be estimated directly. However, in view of 
the fact that these constants have been derived on the basis 
of.a highly idealized structure and loading condition they 
must Бе used with caution. The ratio of the root bending 
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713—717 


1 
"1 + Margin of safety 


_, Ordinate of most highly stressed element 
ii: One-half of wing thickness 


E Actual skin thickness of most highly stressed element 
76 7 Equivalent skin thickness of most highly stressed element 


_ Cross-sectional area òf (assumed) single torsion cell 
Chord X Wing thickness 


_ Perimeter of torsion cell 
Twice the chord 


__ Effective perimeter 7 of torsion cell 
Actual perimeter of torsion cell 


__ Width of equivalent sheet 
i ord 


Average ordinate of upper skin 
^ Maximum ordinate of upper skin 


M Equivalent thickness of lower skin 
m Equivalent thickness of upper skin 


Maximum ordinate of lower skin 
One-half of wing thickness 


"- Average ordinate of lower skin 
1477 


Maximum ordinate of lower skin 


711373 


715— mol nmi? тати?) 


° Allowable torsion stress 
T^ Allowable bending stress 
ШІК; 
UU S cos А 


b\2 
7185-7011 (5) 


7375 
mmis 


Equivalent perimeter 7 of torsion cell 
Actual perimeter of torsion cell 


stiffness to the root torsional stiffness can be estimated 
by means of equations given in appendix B or, preferably, 
from experience with structures similar to that under 
consideration. 

The spanwise stiffness distributions need not be known in 
detail in order to use the charts and approximate formulas. 
If the wing is solid or nearly solid or if its cross sections are 
geometrically similar at all points, the charts for stiffness 
distributions proportional to the fourth power of the chord 
are used. If the wing does not have large cutouts and is 
designed for a constant stress level, the charts for the stiff- 
ness distributions associated with constant stress are used. 
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The use of these charts tends to overestimate aeroelastic 
effects to some extent because, although actual wings are 
designed for constant stress over most of the span, the por- 
tions near their tips are designed on the basis of other con- 
siderations, such as handling loads or minimum standard 
sheet thicknesses; therefore, wings tend to be stiffer near the 
tip than they would be if designed on the basis of constant 
stress throughout. This difference in stiffness is particu- 
Jarly large if the taper ratio is zero. 

If the wing contains large cutouts or if, for any other 
reason, the wing stiffness distribution із known to be sub- 
stantially different from а constant-stress type, the charts 
can be used to furnish semiqualitative results for the various 
aeroelastic phenomena by ‘using fictitious stiffnesses, рго- 
vided the actual stiffness distribution is known at least 
approximately. 'The root stiffnesses of these fictitious dis- 
tributions may be assumed to be the ones that give rise to 
twist or bending angles at the tip which are the same as 
those of the actual wing if the bending moments or torques 
vary as the square of the distance from the tip. For con- 
venience, the spanwise distribution of these fictitious stiff- 
nesses may then be assumed to be proportional to the fourth 
power of the chord. On the basis of these assumptions, 


l _3) -L ads 
ET =з | zr (1—5*)%4в (10) 
where the subscript e refers to Ше fictitious stiffness, and 
where the integral represents the moment of inertia of the 


aren under the function ay plotted against s* about the 


point s*—1. The fictitious torsional stiffness at the root 
can be obtained in the same manner.. The aeroelastic 
phenomena can then be estimated by use of these fictitious 
root stiffnesses and the charts for. stiffness. distributions 
proportional to the fourth power of the chord. 

In the derivation of the charts the wing is considered to 
be clamped at the effective root perpendicular to the elastic 
axis. From the data and analyses presented in references 
1, 2, 4, and 5 & satisfactory location of the effective root 
appears to be at the intersection of the elastic axis and the 
side of the fuselage. 

Tf the rotations at this effective root are known as a 


result of deflection tests or a detailed analysis such ав that. 


of reference 5, the root twist due to torque and the root 
bending due to bending moment may be taken into account 
by moving the effective root inboard by the distance 


две 2 (GJ) (11). 
-% (бл, 
ог 
qu 


where o,7 is the angle of twist at the root due to а root 
torque Т,, and where Г“ is the deflection slope at the 
effective root due to а bending moment М,. Since the 
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distances Ав” and As" may differ from each other, some 
compromise between the two must be made; for unswept 
wings the use of As* appears to be indicated, whereas for 
highly swept wings the use of As" is more appropriate. 


PRELIMINARY SURVEY OF ABROELASTIC BEHAVIOR 


The information contained in some of the subsequent 
sections has been summarized in figure 2 for the purpose of 
ascertaining in advance of more detailed estimates, if desired, 
whether the aeroelastic phenomena considered herein are 
likely to affect the design of the wing structure. This pre- 
liminary survey is not essential to any of the further calcu- 
lations but may show them to be unnecessary in some cases. 

Тһе charts of figures 2 (а) to 2 (d) pertain to wings of 
taper ratios 0, 0.2, 0.5, and 1.0 and constitute plots of the 
dynamic-pressure parameter 4% defined by equation (7) 
against the sweep parameter k defined by equation (9). 
These parameters contain the root stiffnesses (GJ), and 
(ЕТ),; 1, when а preliminary survey of aeroelastic effects is 
to be made, no information whatever concerning the wing 
stiffness is available, the following relations for 4% and k 
may be used: 

Сі, E (1--X)* cos ^ | 
(13) 


в сапы 1! 
18432 g g, пай, 
т 1+0 9 A, 
k= LEA Е quei ль tan A (14) 


where F. is а root-stiffness parameter defined in appendix 
В and plotted in figure 3, and where та and т are defined by 


__ 13767617" 


(15a) 
7417377971157 17 
апа 
ызды Мез (15b) 
7157187197] 16 


in terms of some of the factors defined in table 1. 

Figure 2 (e)Jpertains to wings for which the moment arm 
e is zero and, hence, k is infinite; with the degree of approx- 
imation involved in the use of figures 2 (а) to 2 (d), figure 
2 (e) can be used for wings with |k] > 25. This figure con- 
sists in а plot of the dynamic-pressure parameter 9, defined 
by equation (8), against the taper ratio Х. If no informa- 
tion is available concerning the root bending stiffness (ZI), 
contained in 4, the following relation may be used: 


a (l HAY Ах sin А 
(16) 


Тһе various lines of Ше charts of figure 2 designate the 
conditions at which а wing designed on Ше basis of strength 
considerations alone is likely to encounter divergence or span- 
wise shifts of the center of-pressure by various amounts; 
positive shifts are those toward the tip. These spanwise 
shifts furnish an estimate of the increase in root bending 
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(а) Wingg of taper ratio 0. 
FIGURE 2.— Charts for a preliminary survey of aeroelastic phenomena. 


phenomena of unswept wings are similar to those of swept- 
The &erodynamic-center shift asso- 
ciated with the shift in the lateral center of pressure у or in 


forward wings if e; is positive апа to those of sweptback 


wings if e; is negative. 


except 


for small positive values of k (associated with-sweep angles 
А negative value of e; 


may be obtained at supersonic speeds, but under most con- 
Similarly, jn figure 2 (e) 


a negative value of 4 implies that A is negative (that is, that 
whereas 8, positive value implies 


The significance of negative values ‘of g* is that e; is nega- 
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In using figures 2 (а) to 2 (d), estimates must be made of 
either the root stiffnesses (jn conjunction with eqs. (7) апа 


that the wing is swept back. 


The factor F, is obtained from figure 3 for 


(9)) or of the effectiveness factors y, and m (for use in eqs. 
the largest value of e; likely to be encountered at the design 


(13) and (14)). 


q* is calculated for the combination of dynamic pressure q, 
lift-curve slope Cz, and moment arm e, which is likely to 


load factor and for the given taper ratio А. The parameter 
be critical from an aeroelastic point of view. For an unswept 
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certain sweep 


In other words, the aeroelastic 


For unswept wings k is approximately equal to zero, and Ше 


üeroelastic phenomena referred to in the charts of figures 
2 (а) to 2 (d) are similar to those of swept wings defined by 


points in quadrant 2 if q* is positive and by points in quad- 


rant 3 if 4% is negative. 
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FiaunE 2.—Continued. 


(b) Wings of taper ratio 0.2. 


aeroelastic effects on the spanwise center of pressure or the 


possibility of an approach to the divergence condition, further 


The charts of this report may be 
used for the preliminary calculations; once Ше structure has 


-calculations are desirable. 


proceeds in the 


then 
same manner ав for figures 2 (в) to 2 (d) using, in the parameter 


a * 
в 


If the moment arm е, is so small or Ше angle of sweep so 
large that the parameter Ё exceeds the range covered by 
figures 2 (a) to 2 (d), the chart of figure 2 (e) may Бо used for 
the purpose of a preliminary aeroelastic appraisal of the given 
wing. Іп this figure only the parameter 9 is required, since 
Е is considered to be infinite. The parameter 9 can be ob- 

the condition for maximum оС, 


been designed, more refined methods such as that of reference 


1 may be used. 
tained from equation (16). The analysis 


4, 


wing Ше combination for which the product 4Сг, бі is а 


is & maximum is likely to be critical. 


а, 


maximum is likely to be critical; for a swept wing the combi- 


nation for which gC; 


On the other hand, 


If the shift in span wise center of pressure (and 
if the point on the chart indicates the likelihood of significant 


The values of g* and Е define a point on one of the charts 
of figures 2 (а) to 2 (d) (whichever is closest to the actual 


taper ratio). 
any associated shift in the aerodynamic center) at that point 
is small and, in the case of ап unswept ог a sweptforward 


wing, if the absolute value of the ratio of the value of q* at 
that point to the value of 4% at divergence for the value of 


k at that point is small, the static aeroelastic phenomena 


The parameter k is then calculated for the same value of е. 
discussed in this report probably need not be taken into 


account in designing the wing structure. 
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(c) Wings of taper ratio 0.5. 
Етоонв 2.—Continued. 
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(d) - Wings of taper ratio 1.0. 
FIGURE 2.— Continued. 
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Етаовв 2.—Concluded. 
Етсовя 3.— The root-stiffness function Fn 
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CALCULATION OF THE VARIOUS AEROELASTIC PHENOMENA 


Dynamic pressure at divergence.— Тһе solutions for the 
divergence speed obtained by the direct method in reference 
2 and those obtained by the numerical matrix method given 
in appendix А сап be summarized by approximate formulas 
which give the dimensionless parameters g*p and gp (the 
values of the parameters defined in eqs. (7) and (8) that 
correspond to the value of the dynamic pressure g at diver- 
gence) in terms of their ratio Ё defined by equation (9). 

These approximate formulas are 


“РК (18) 
and 
nS K 
dp-————— (19) 
K:— 


Wher the angle of sweep is zero, equation (18) reduces to 
Фе Ка, and when the moment arm e, is zero, as it may be 
in supersonic flow, equation (19) reduces tò 4р-- E The 
2 
constants А, and Ко are given in table 2 for wings with taper 
ratios of 0, 0.2, 0.5, and 1.0 for both types of stiffness distribu- 
tions; the parameter q*» for unswept wings and the parameter 


Gp for swept wings with e;—0 are plotted against the taper. 


ratio А in figures 4 (a) and 4 (b), respectively. 

With the values of 4“ or Gp given by equations (18) and 
(19) and the definitions of these two parameters given by 
equations (7) and (8), the values of 4 required for divergence 
may be determined. If desired, the corresponding airspeed 
may be determined from the relation 


= / QD 


The value of gp is often negative for sweptback wings, and, 
since a negative dynamic pressure does not correspond to any 
real speed, these wings cannot diverge. These negative 
values of gp, nonetheless, are useful as reference values in 


TABLE 2.—VALUES OF THE COEFFICIENTS К, AND К, 


К K3 
1 


Ву ana- Ву апа- 
Бах ytlo inte- DEIN уйе Inte- 
А оп п ntegration on 
fret 2) fret. 2) 
2.58 Ё 
2.83 5 А 


2.47 
2.74 480 
2.81 
2.25 
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other aeroelastic phenomena. (If the negative value of до 
is very low in absolute magnitude, divergence in a higher 


mode associated with a positive qp is conceivable. However, 
3.07 „Stiffness proportional to c^ (matrix integration) 
Иа 
Stiffness proportional to сё (analytic integration}- oe 
| 20 | Stiffness given by constant-stress criterion for (ИУ - |.0 
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О 5 1.0 15 
Taper ratio, А 
(а) Effect of taper ratio on 4*р for unswept wings. (g*p-- Кі in this 
caso.) 


(b) Effect of taper ratio on gp for е1-+0. (>= -5 in this case.) 


Етсове 4.—Effect of taper ratio on the dynamic pressure at divergence. 
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for the ordinary wings with straight leading and trailing 
edges and substantially straight elastic ахев to which the 


analysis of this report is applicable, the magnitude of the 


lowest positive до is always much larger than the absolute 
magnitude of the negative qp associated with the first mode. 
Higher-mode divergence of а practical wing for which а neg- 
ative до is indicated by the analysis of this report is therefore 
very unlikely.) 

The values of the constants K, and К, given in — 2 
differ somewhat from the corresponding values resulting from 
the matrix solution in appendix A. The matrix results are 
probably more significant because they are based upon more 
realistic aerodynamic assumptions; the А! and К» values in 
reference 2 tend to give conservative results. 

The value of др calculated for any given value of 4%р or 
qp depends on the value of the effective lift-curve slope C. 


ог cı, and, hence, on the Mach number. Аз suggested in 


references 1 and 2, the value of др calculated at various 
Mach numbers may be plotted against Mach number. If 
lines of the actual dynamic pressure at several altitudes as a 
function of Mach number are drawn on the same plot, an 
intersection of the divergence line with one of the lines of 
actual dynamic pressure designates possible divergence at 
that value of dynamic pressure, Mach number, and altitude. 
If this plot is on log-log coordinates, the lines of actual dy- 
namic pressure are straight and the ratio of the dynamic 
pressure at divergence to the actual dynamic pressure at а, 
given Mach number and altitude can be scaled off directly. 
(See refs. 1 and 2.) 

Spanwise angle-of-attack distributions.—In appendix А, 
an approximate expression is determined for the change in 
angle ‘of attack due to wing flexibility. For the additional- 
type angle-of-attack distribution (a, is constant) the angle 
of attack due to structural deformation a, is given by - 


а 1 g/g» Е, А 20 
көгім fy) | (20) 


ор R 


4р 


The functions у, and Af, of the spanwise coordinate s* and 
the function F, of the parameter Г are given in figure 5 for 
swept wings with taper ratios of 0.2, 0.5, and 1.0 and with 
the two different types of stiffness distributions. For wings 
with zero taper ratio the structural deformation cannot be 
obtained from equation (20), as is pomted out in appendix 


A. However, the ratio к > аз а function of the spanwise 


coordinate 8* is shown in figure 6 for the two different stiff- 
ness distributions, several values of 4/0р, апа several values 
of the parameter k. 

The spanwise distribution of o, due to а linear twist 
(аҙ--8%о,), which may be either symmetric or antisymmetric, 
13 approximately 


90 (6,4 p Af; (21) 
p [ws 


до 


Og _ 


757 


where the functions №, Аў, and Р, are given in figure 7 for 
wings of taper ratios 0.2, 0.5, and 1.0. Тһе angle-of-attack 
ratio is shown in figure 8 as a function of в*, 4/4р, and Е for 
wings of zero taper ratio. 

The results of equations (20) and (21) may Бе superim- 
posed. For example, if the spanwise distribution of a, due 
to rolling is to be found, these equations must be added in 
such proportion that 


ож pb 
Put 
EN 
Ур tF 
so that 
b, b »\ 26 | 


where pb/2V is the angle of attack аф the tip due to roll. 
The distribution of o, due to roll 1s then 


^ = (1— S) (НРА - Gs FAT) | (23) 


2V Ор 


Spanwise lift distributions.—If desired, the lift distribu- 
tions can be obtained for the angle-of-attack distributions 
given in the preceding section by one of the commonly used 
methods of calculating spanwise lift distributions, such as 
that of reference 6 or by the approximate method of reference 
7. However, the following method is simpler and, within 
the approximation of the present report, as accurate рго- 
vided the rigid-wing (4--0) loading J, is calculated by Ще 
methods of reference 6 or 7 or is obtained from the charts of 
reference 8 or, in general, by an accurate analytical method. 

Within the framework of the assumptions made in the 
analysis the lift per inch of span is proportional to the local 
angle of attack, so that 


potted (24а) 


for geometrical angles of attack.which’are constant along 
the span, and 


n ------ (24b) 


for geometrical angles of &ttack due to linear twist, where 
55 and о/о, are obtained as indicated in the preceding 
section. 

If no better approximation is available for the loading lo 
(which is likely to be the case at transonic speeds), it may 
be estimated for geometric angles of attack which are con- 
stant along the span from the relation 


іе: Оһ, саа, (25а) 


and for all other geometric angles of attack from the relation 


EUNTES (25b) 
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FIGURE 5.— The angle-of-attack distribution functions Л, АЛ, and F, for constant geometric angles of attack. 


(а) Stiffnesses given by constant-stress criterion for 
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(b) Stiffnesses related to those given by constant-stress criterion for wings with taper ratio 0.5, 


. Fraunz 5.— Continued. 
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(c) Stiffnesses proportional to сї. 
FraunE 5.—Concluded. 
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(а) Stiffnesses proportional to сі. | (b) Stiffnesses given by constant-stress criterion. 
Етатнк 6.— The angle-of-attack distributions for constant geometric angles of attack for wings of taper ratio 0. 


REPORT 1140— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


762 


x 


j Q | 
LNALLLLLISAL IT LL: 


А6 
ОЗГЕ АГК 


"m | 
XL 


clion A fox 


Function № 
" 
— = 


ККЕ ЕГА 


с) © © © “ с © q 


Фу рио ©, ' ; 


о 


о 


(hic) 4. 
(Ме); и 


Етатвв 7.— The angle-of-attack distribution functions fs, Аў», and F3 for linearly varying geometric angles of attack. 


(а) Stiffnesses given by constant-stress criterion for 
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(b) Stiffnesses related to those given Бу constant-stress criterion for wings with taper ratio 0.5. 
Етаове 7.—Continued. 
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(c) Stiffnesses proportional to с“. 
Етсовв 7.—Concluded: 
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(b) Stiffnesses given by constant-stress criterion. 


(а) Stiffnesses proportional to ct. 


FraunE 8.—The angle-of-attack distributions for linearly varying angles of attack for wings of taper ratio 0. 
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Lift and moment coefficients—The wing lift coefficient 
Съ the wing-root bending-moment coefficient Св, and the 
wing-root twisting-moment coefficient Ст шау be obtained 
in terms of their respective rigid-wing values by means of 
the following approximate expressions: 


now (1— 
Cr, 1 "ue y) 


= (26) 
бла 19 
до 
4 
1——(1—py) 
оь 
О», ML t 
до 
an (1 — 
б... 1 25 (1--т») (28) 
Ст, ee A 
dp 


where the coefficients », и, and + depend on the type of 
loading. The subscript.1 is used for additional-type angle- 
of-attack distributions and the subscript 2, for linear-twist- 
type angle-of-attack distributions. The coefficients и, ш, 
and ті are given in figure 9 as functions of the parameter k 
for wings of taper ratios 0.2, 0.5, and 1.0. For wings with 
taper ratio zero the ratios of the lift, bending-moment, and 
twisting-moment coefficients to their respective rigid-wing 
values are given in figure 10 as a function of 4/4 for several 
values of the parameter k. The values of уз, ш, and та are 
given in figure 11 for wings of taper ratios 0.2, 0.5, and 1.0, 
and ratios of the lift, bending-moment, and twisting-moment 
coefficients are given in figure 12 for wings of zero taper ratio. 
The additional-twist and linear-twist results of equations 
(26) to (28) may be superimposed in the same way as those 
of equations (20) and (21). 

The wing rolling-moment coefficient C;, is defined as the 
rolling moment of the loads on both wings about the fuselage 
center line divided by gSb. Therefore, 


2M, cos A--2T, sin A+5 То 
0. a8 89 
The angle-of-attack distribution due to rollmg given in 
equation (22) must be used in finding the values of М,, 7;, 
and Ly in equation (29). | 
Spanwise centers of pressure and aerodynamic centers.— 
The spanwise location of the center of pressure is given by 
the distance 


o0 Св 
8—5 E (30) 
or the dimensionless distance 
ЕРКІ 
“oC (81) 
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1 
(Inasmuch as 7* is equal to 8“ by virtue of the definitions of 
those dimensionless quantities (see also fig. 1), eq. (81) can 
be considered to be an expression for 7* rather than 8%, if 
desired.) With the values of the bending-moment and lift 
coefficients given in the preceding section, the ratio of 8 to 
its value for the rigid wing may be calculated from either of 
the equations 


і 


X tn 
аса 
80 4 
les (b 
| 2-01) 
(82) 
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where the shift in spanwise center of pressure А5 is defined 
as 5— 50, and where ш and >; are used for constant geometrical 
angles of attack and p: and уз, for linearly varying geometri- 
cal angles of attack. 

The shift due to aeroelastic action of the longitudinal 
position of the center of pressure associated with a given 
shift of the spanwise center of pressure is equal to sin A 43. 
The shift in aerodynamic center (positive when rearward, or 
stabilizing) can consequently be calculated by substituting 
into equation (17) the values of Аз obtained from equation 
(32) with values of ш and и. 

Inertia effects.—No charts are presented in this report for 
the effects of inertia on quasi-static aeroelastic phenomena— 
that is, aeroelastic phenomena associated with flight at 
constant acceleration; the manner in which mass is distrib- 
uted varies so widely among different wings that preparation 
of a generally applicable set of charts for inertia effects 
appears to be impractical at present. Furthermore, except 
for flying wings, the wing deformations due to inertia loads 
are small compared with those due to aerodynamic loads, 
the two types of loads being in about the same ratio as the 
wing weight to the weight of the entire airplane. If desired, 
however, inertia effects and the aeroelastic increment in 
these effects can be estimated in the manner described in 
the following paragraphs. 

From the known or estimated mass distribution of the 
wing the inertia load J, per inch of span and the inertia torque 
tı per inch of span can be calculated for any given normal, ` 
pitching, or rolling acceleration. Substitution of these loads 
and torques for the terms | and Фес in equations (АЗ) ог 
(A36) and equations (А2) or (A35), respectively, yields the 
values of the accumulated bending moment and torque due 
to the distributed inertia loads and torques. Іп turn, sub- 
stitution of these accumulated bending moments and torques 
in equations (А4), (Аб), and (Аб), or in equations (А37) 
and (А38) and the matrix equivalent of equation (Аб), 
yields the angle-of-attack distribution due to the deforma- 
tions caused by the inertia effects associated with the given 
acceleration. 
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(c) Stiffnesses proportional to с“. 


ЕтаовЕ® 9.—Concluded. 
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(b) Stiffnesses given by constant-stress criterion. 


Frayre 10.— The lift- and moment-coefficient ratios for constant geometric angles of attack for wings of taper ratio 0. 


(а) Stiffnesses proportional to сі. 
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Fiaunzs 11.—The lift- and moment-coefficient parameters по, го, and та for linearly varying geometric angles of attack. 
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(b) Stiffnesses given by constant-stress oriterion. 


Fraunx 12.— The lift- and moment-coefficient ratios for linearly varying geometric angles of attack for wings of taper ratio 0. 


(а) Stiffnesses proportional to c*. 
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Тыв angle-of-attack distribution can be considered as а 
geometrical angle-of-attack distribution. For the purpose of 
calculating the increment caused by aeroelastic action, this 
distribution сап be approximated by а linear-twist angle- 
of-attack distribution with & value at the wing tip which 
18 such that the moment about the effective wing root of 
the area under the linear-twist distribution equals the mo- 
ment of the area under the calculated angle-of-attack dis- 
tribution due to inertia effects. (The moment, rather than 
the area, is suggested as a basis of correlation because the 
angles of attack near the wing tip are more important. in 
aeroelastic phenomena than those at the wing root.) The 
justification for this rather arbitrary approximation to the 
angle-of-attack distribution 18 as follows: As previously 
mentioned, the wing deformations due to inertia loads are 
likely to be small compared with those due to aerodynamic 
loads; furthermore, the correction to be applied to these 


deformations as a result of aeroelastic action is usually | 


small compared with these deformations and, hence, is very 
small in comparison with the total wing load, so that the 
correction need not be calculated as accurately as the cor- 


rection for aeroelastic effects to the rigid-wing lift distribution. - 


The angle of attack due to structural deformation a, asso- 
ciated with the linear-twist distribution can then be obtained 
from. equation (21) and figure 7 or, if ^=0, from figure 8. 
The lift distribution associated with the total angle-of-attack 
distribution due to the deformations caused by the inertia 
effects, including the increment in this angle-of-attack dis- 
tribution produced by aeroelastic action, can then be found 
from equation (24b), in which a, and l pertain to the cal- 
culated angle-of-attack distribution due to the inertia effects 
(not the linear approximation to this distribution). This lift 
distribution can be integrated to obtain the lift, bending 
moment, rolling moment, and aerodynamic-center position 
due to inertia effects, as modified by aeroelastic action. 

The lift and rolling moment calculated in this manner 
may then be combined with the lift and rolling moment for 
steady. level or rolling flight calculated by the method out- 
lined in the preceding sections. For instance, if the con- 
tributions of the tail and the fuselage to the airplane lift 
can be neglected, the wing lift can be written as 


п) Cra ets EE) п 


where (5 = 
n 
due to inertia effects, including aeroelastic effects; it is equal 
to —W, plus the lift on both wings due to inertia effects, 
as modified by aeroelastic action, per unit load factor and 
is almost always negative. In the preceding equation Cha, 
is а wing lift-curve slope which includes static aeroelastic 
effects and is equal to Cz, multiplied by the factor on the 
right side of equation (26). Then | 


1 А 1 Ор, 245 
"144, т (8 o] W—W, 
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1 Che, 249 


is the total normal force per unit load factor 
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where 
1 


Оз, i 1 Cha, 


je ll A) 
W-— ТР» on 2 


is а wing lift-curve slope which includes static aeroelastic 
effects, inertia effects, and aeroelastic modification of the 
inertia effects. 


ILLUSTRATIVE EXAMPLE 


The parameters of a swept wing, which differs from the 
wing of the illustrative example of.reference 1 only in the 
width of the fuselage to which it is attached, are given in 
table 3. The values of As? and As" were calculated from 
the dimensionless root-rotation constants used in the 
example of reference 1, Фу <0 ара Qr,,— —0.25, by means 
of the relations 

As*= Qo, We 


Ав == Qr We 


1 


where w,, as defined in reference 1,is the distance along the 
span between the effective root and the innermost complete 
section of the torsion box perpendicular to the elastic axis. 
In the wing of the illustrative examples of the present 


TABLE 3.—PARAMETERS OF WING USED IN 
ILLUSTRATIVE EXAMPLE 


Geometrical parameters 
A : ада 


былды абы ыты шаш --0.444 
(87) 16а іп, ...-.......8.04Х10 


CET», Ibjsq іш. 7. 9.58% 


-——— ————— —. ——————.—.. 
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report and of reference 1, w, із 22.4 inches. These relations 
for As? and As" can be obtained from equations (11) and (12) 
of the present report in conjunction with the definitions of 
the root-rotation constants given in equations (158). and 
(15d) of reference 1; in the notation of the present report 
the definitions are 


— er l1, 
Geo GJ); 
Qr = T,MIM, 

м VJ (E), 


The stiffness 15 assumed to vary as the fourth power of the 
chord in the example of the present report. 

The subsonic and supersonic values of the parameter Ё 
were calculated from equation (9). By means of ap- 
propriate values of the constants A, and К» taken from table 
3, the values of g*p were calculated from equation (18) and 
included in table 3. From these values of g*5, Ше sub- 
sonic and supersonic dynamic pressures at divergence were 
found through the use of equation (7) and are given in table 
3. These values of gp vary as the reciprocal of the effective 
lift-curve slope, the corresponding values of e; being assumed 
to remain constant. 

In order to find the angle-of-attack distribution for 
additional-type loadings from equation (20), the values of 
Е, and of the functions f and Af, were taken from figure 
5 (с). The spanwise change in angle of attack is shown in 
figure 13 for different values of the dynamic-pressure ratio. 
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FIGURE 13.— Effect of aeroelastic action on some aerodynamic proper- 
ties of the wing used in illustrative example. 
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The values of rı, ті, and и, were obtained from figure 9 (с) 
and substituted into equations (26), (27), and (28). The 
wing lift coefficient, wing rolling-moment coefficient, and 
spanwise center-of-pressure ratios, as well as the shift in 
aerodynamic center, were calculated by use of these ap- 
proximate equations in conjunction with equations (17) 
and (29) and are shown in figure 13 as functions of the 
4 
4 


dynamic-pressure ratio — 


г DISCUSSION 


1 


LIMITATIONS OF THE CHARTS AND APPROXIMATE FORMULAS 


The charts and the approximate formulas presented in 
this report are subject to certain limitations as a result of 
the approximations made in the calculations on which they 


 &re based. These limitations take the form of restrictions 


on the plan form, on the speed regime, and on the wing 
structure. The results obtainable by the use of the charts 
are likely to be unsatisfactory for wings of very low aspect 
ratio or very large sweep and relatively unsatisfactory for 
wings of zero taper ratio. 

Wings of low aspect ratio are ruled out on three counts: 
(1) the extent to which aerodynamic forces are overestimated 
in replacing the wing by one with an effective root and tip 
is larger for wings of low aspect ratio than for wings of 
high aspect ratio, (2) elementary beam theory is unsatisfac- 
tory for calculating the deformations of wings of very low 
aspect ratia (because the effects of end constraint, shear lag, 
shear deformation, and bending-torsion interaction are 
more important when the aspect ratio is low), and (8) the 
assumptions made concerning the lift distribution of the 
Wing are more nearly true for wings of high than for those 
of low aspect ratio. | 

For wings with very large angles of sweep, also, the use of 
an effective root and tip introduces relatively large errors in 
the aerodynamic forces. Furthermore, the root rotations 
neglected in the calculations (bending rotation due to tor- 
sion and twist due to bending) are likely to be important 
for wings with large angles of sweep. 

The aeroelastic analysis of wings with zero taper ratio 
entails certain mathematical difficulties which do not arise in 
the case of wings with nonzero taper ratio. The stiffness of 
such wings is zero at the tip and very low near the tip, so that 
the boundary conditions for о and Г given by equations 
(А10а) to (А10с) in appendix A are indeterminate. Аз а 
result of the relatively large values of the reciprocal of the 
stiffness near the tip, the numerical-integration methods used 
in the matrix calculations are less accurate. These difficulties 
also occur in other methods of solving the aeroelastic equa- 
tions, such as energy methods. Furthermore, the structural 
behavior near the wing tip is not represented adequately by 
elementary beam theory. Finally, that the aeroelastic results 
calculated for wings of zero taper ratio are not as reliable as 
those for other wings is evidenced also by the fact that they 
do not lend themselves to systematization by means of 
approximate formulas, as do the aeroelastic results calculated 
for other wings. 

As a result of these considerations delta wings are unsuit- 
able for aeroelastic analysis by means of these charts because 
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they have low aspect ratios, large angles of sweep, and zero 
taper ratio. 

In order to use the charts two aerodynamic parameters 
must be known for any given case: the effective wing lift- 
curve slope and the section aerodynamic center. From an 
aerodynamic point of view the charts of this report may be 
used in almost all cases for which thesé quantities are known. 
The exceptions stem from the fact that the spanwise distribu- 
tion of the lift on the rigid wing is assumed to be proportional 
to the chord, and the distance from the section aerodynamic 
center to the elastic axis (as a fraction of the chord) 18 
assumed to be constant along the span. These assumptions 
are not valid for wings with large angles of sweep and wings of 
low aspect ratio, аз implied previously. They are also invalid 
to a greater or lesser extent for most wings in the transonic 
region. Consequently, even when the lift-curve slope and 
the section aerodynamic center are known, any results 
calculated for transonic speeds must be used with caution. 
» Another aerodynamic assumption implied in the charts is 
that no concentrated aerodynamic forces, such as those due 
to a tip tank or nacelle, act оп the wing. Relatively small 
nacelles in the inboard half of the span can probably be 
ignored for the purpose of an aeroelastic analysis at the 
preliminary design stage. However, large tip tanks cannot 
usually be ignored even in a preliminary aeroelastic analysis; 
the aeroelastic phenomena may in such cases be greatly 
underestimated by calculations made with the charts of this 
report. 

The assumption concerning the applicability of elementary 
beam theory to the calculation of wing deformations due to 
neroelastic action serves to restrict the wings that can be 
analyzed by means of the charts to those of moderate or high 
aspect ratio, as stated previously. Neglect of chordwise 
bending (elastic camber) effects in the calculations on which 
the charts are based serves to:impose a lower limit on the 
thickness of the wings for which the charts may be used. 
Whether this limit is within the region of practical thicknesses 
is questionable, however. The divergence tests of reference 
2, which were performed on flat plates of moderately high 
aspect ratio and with a thickness of 2.5 percent, showed no 
obvious chordwise bending effects, although the relatively 
small differences between the measured and calculated 
divergence speeds may have been due in part to such effects. 

As mentioned previously, for wings with taper ratios 
between 0 and 0.2 the results of aeroelastic calculations are 
likely. to be relatively unreliable. For taper ratios greater 
than 0.2, the stiffness of actual wings tends to be greater 
near the tip than that given by the constant-stress criterion; 
consequently, any given aeroelastic effect is likely to be some- 
what less than that calculated on the basis of а constant- 
stress stiffness distribution but much larger than that 
calculated on the basis of a c* distribution. 


If а given structure contains large cutouts which give 


rise to discontinuities in the stiffness distributions, equation 
(10) can be used to calculate a fictitious root stiffness to be 
used in conjunction with charts for c*-type stiffness distri- 


butions, provided the magnitudes of the discontinuities аго 


known or can be estimated. 
Use of the charts of this report is premised on the assump- 
tion that the elastic axis is at ап approximately constant 
321000 —00———49 
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fraction of the chord. If the location of the elastic axis 
varies somewhat along the span, the use of an average value 
tends to give satisfactory results for the aeroelastic phe- 
nomena of swept wings; for unswept wings, however, the 
results obtamed on the basis of this approximation have to 
be used with caution. If the elastic axis exhibits abrupt 
shifts along the span as a result of large cutouts or for other 
reasons, the charts should be used only for moderately or 
highly swept wings. This restriction is mitigated to a 
certain extent by the fact that an abrupt shift in the locus 
of shear centers does not necessarily imply an equally large 
or equally abrupt shift in the elastic axis. 


RELATION BETWEEN STRENGTH AND STIFFNESS AS DESIGN CRITERIA. 


The strength of a structure is its ability to withstand 
applied loads without failure; the stiffness of a structure is 
its ability to deform relatively little under the applied loads. 
The two terms are related (a fact which forms the basis of 
the constant-stress type of stiffness distributions used in 
this report) but are not synonymous. The problem of when 
to design for strength and when to design for stiffness and 
the related problem of how to design a wing for stiffness 
when required to do so have been recognized for a long time. 
Because of the complexity of these problems no generally 
satisfactory solution exists at present, but the charts pre- 
sented herein shed a certain amount of light on the problem 
insofar as stiffness requirements occasioned by the aero- 
elastic phenomena considered in this report are concerned. 

The charts of figure 2 indicate the extent to which a wing 
is likely to be affected by aeroelastic phenomena—that is, 
how far it is from divergence and how much its spanwise 
center of pressure is likely to shift as a result of aeroelastic 
action, provided the wing is designed on the basis of strength 
considerations alone. Jf the margin against divergence is 
too small, or if the spanwise center of pressure and the 
associated shift in the aerodynamic center are deemed 
excessive, the wing has to be stiffened beyond the amount 
associated with the required strength. The charts of figure 2 
therefore serve to delimit the regions in which a wing can 
be designed on the basis of strength considerations alone 
and those in which stiffness considerations predominate, 
at least to the extent of satisfying the stiffness requirements 
associated with the aeroelastic phenomena considered herein. 

The bending moment of inertia required by considerations 
of strength alone for the root section of a wing is directly 
proportional to the design load n(W—W,,), to the spanwise 
coordinate of the center of pressure, and to the wing thick- 
ness at the root and is inversely proportional to the allowable 
bending stress Ер. Alternatively, this bending moment 
of inertia may be considered to be proportional to the design 


wing loading ue e) Lin We) 


the wing thickness ratio at the root, and: to a function of 
14-2Х 


› to the square of the wing area, to 


_ the taper ratio (which is + ИА if strip theory is assumed to 


apply) ;itis inversely proportional to Ев and independent 


of the aspect ratio. These relations for the bending 
moment of inertia required by considerations: of strength 
alone аге implied in the following discussion of structural 
requirements imposed by considerations of stiffness, 


N 
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In general, a wing with а high value of 4% (see eqs. (7) and 
(13)) is most likely to be affected by aeroelasticity (вее fig. 2) 
апа, for а given value of 4%, swept wings аге much more 
likely to be affected by aeroelasticity than unswept ones. 
(See fig. 2 and eqs. (9) and (14).) Consequently, the follow- 
ing wings are most likely to be subject to aeroelastic phe- 


nomena, provided they are designed | on the basis of strength 


considerations alone: 

(1) Wings operating at a high flying speed or high ner 
pressure 

(2) Swept wings 

(3) Thin wings 

(4) Wings designed for a low wing loading 

(5) Unswept and moderately swept wings with an elastic 
axis relatively far back on the chord or likely to fly in a 
condition in which the section aerodynamic centers are 
relatively far forward on the chord 

(6) Wings operating at a Mach number at which the lift- 
curve slope is relatively high 

For given wing loadings and given wing areas, some 
aeroelastic phenomena of wings designed on the basis of 
‚ strength considerations alone appear to be substantially 
unaffected by changes in the taper ratio—for imstance, the 
spanwise shift of the center of pressure and the dynamic 
pressure required for divergence. (In the case of the dynamic 
pressure required for divergence, the parameter q*p (fig. 4), 
the root stiffness, and the root chord decrease with increasing 
taper ratio, and the net effect of taper is small.) Ор the other 
hand, the change in the lift due to aeroelastic action is more 
sensitive to the taper ratio; it is more significant for wings 
with high taper ratio than for wings with low taper ratio. 

The effect of aspect ratio on aeroelastic phenomena tends 
to be small for unswept wings of a given wing area, because 
these phenomena are determined largely by the magnitude of 
the parameter 4%, which is independent of the aspect ratio for 
а, given wing area. Жог the aeroelastic phenomena of highly 
swept wings, however, the parameter 4 is more significant. 
This parameter is proportional to the swept-span aspect 
ratio for wings of a given area. Consequently, with a given 
wing area, taper ratio, and design wing loading, the аего- 
elastic effects of swept wings tend to be more pronounced for 
wings with high aspect ratio than for those with low aspect 
ratio. This statement is particularly true for the shift of the 
aerodynamic center, because a given spanwise shift of the 
center of pressure résults in a much greater chordwise shift 
in the case of a swept wing of high aspect ratio than in the 
case of a swept wing with low aspect ratio. 


STRUCTURAL WEIGHT ASSOCIATED WITH THE REQUIRED STIFFNESS 


When а given wing has been shown to be subject to un- 
desirably large aeroelastic effects (by means of the charts of 
this report or by any other method), the problem arises how 
to distribute the additional required stiffness. If, for 
instance, the dynamic pressure on an unswept wing is within 


10 percent of the dynamic pressure required for divergence 
and a margin of 20 percent is desired, an increase of 10 percent 
in the torsional stiffness along the entire span will produce 
the desired result. The question remains, however, whether 
structural weight can be saved by increasing the stiffness 
more than 10 percent in some places and less in others. 

Some insight into this problem may be gained, at loast 
insofar as the aeroelastic phenomena considered herein are 
concerned, from aeroelastic and weight calculations that have 
been made for a family of somewhat arbitrarily selected 
stiffness distributions which differ from the distribution 
required by the constant-stress criterion. The ratios of the 
local stiffnesses to those associated with constant stress are 
shown at the top of figure 14. The structural-weight factor 
Е. is shown for two of these stiffness distributions ав a 
function of the taper ratio. The function Е is proportional 
to the weight W, of the primary load-carrying structure and 
depends on the manner in which the wing stiffness and thick- 
ness are distributed along the span. (See appendix B.) 

The results of the aeroelastic calculations for wings with 
taper ratio 0.5, constant wing thickness ratio h/c along the 
span, and these two stiffness distributions are included in 
table 2 and figures 5 (b), 7 (b), 9 (b), and 11 (b). The 
designation “ехсевв strength" in these figures refers to the 
stiffness distributions increased over the constant-stress 
requirement, as shown in figure 14, with а value of о=2.0. 
The results of the aeroelastic calculations for the stiffness 
distributions decreased below the constant-stress require- 
ment to a value of w=0.5 are the same as those for the 
constant-stress stiffness distributions for wings with linearly 


varying wing thickness ratio and (с Jo, 99. 


The results of the weight calculations and the aeroelastic 
calculations may be combined in several ways. The dynamic 
pressure at divergence, for instance, can be varied by changing 
the bending and torsional stiffnesses uniformly along the 
span, by leaving the stiffnesses at the root unchanged and 
varying the stiffness distribution in a manner similar to that 
indicated at the top of figure 14, or by a combination of the 
processes. А specified dynamic pressure at divergence can 
therefore be obtained as the result of many combinations of 
root stiffnesses and stiffness distributions. Figure 15 (a) 
consists in essence of a plot of the structural weights 1880- 
ciated with combinations of this type against the tip stiffness 
ratio w for a specified dynamic pressure at divergence. This 
figure indicates that the least weight is associated with 
values of the tip stiffness ratio greater than 1. Similarly, 
figures 15 (b) and 15 (с) consist in essence of plots of the 
structural weights associated with various combinations of 
root stiffnesses and stiffness distributions required for shifts 
of +10 percent in the spanwise center of pressure at a 
specified dynamic pressure. Figures 15 (b) and 15 (c) also 


‘indicate that the structural weight is least for values of tho 


stiffness ratio w greater than 1. 
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Тір stiffness ratio, ш 
(а) Relative weight required for а given dynamic pressure аб divergence. 


(b) Relative weight required for 2.0.9 at a given dynamic pressure. 
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(c) Relative weight required for 5—11 аф а given dynamic pressure. 
89 
Fiaure 15.— The effect of the tip stiffness ratio on the structural 


weight required for & given divergence dynamic pressure or given 
shifts in the spanwise center of pressure of wings with taper ratio 0.5. 


The significance of figure 15 is that, if a given wing designed 
on the basis of strength alone needs to be stiffened for 
aeroelastic reasons, most of the stiffening material should be 
added in the outboard regions, provided the weight of the 
material other than that of the primary load-carrying struc- 
ture is unaffected by the stiffening process. Іп fact, on the 
basis of aeroelastic considerations alone, weight might be 
saved in some cases by removing material from the root and 
adding material at the tip; needless to say, however, strength 
requirements would be violated by this procedure. Just 
where the material should be added in'the outboard regions 
cannot be said on the basis of the calculations made for 
figure 15, since these calculations assume any modifications 
to the constant-stress stiffness distributions to be made as 
indicated at the top of figure 14. However, it appears 
unlikely that great weight savings can be had by using modi- 
‘fications which. differ substantially from those of figure 14. 
SOME REMARKS CONCERNING THE AEROISOCLINIC WING 


The term ''aeroisoclinie" refers to wings which deform 
under an aerodynamic load in such а fashion that the angles 
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of attack of all sections relative to the free stream remain 
unchanged. Such a wing has the advantage that its аего- 
dynamic loads do not change under aeroelastic action either 
in magnitude or in distribution; its aerodynamic center, for 
instance, is unchanged, and the wing cannot diverge. Tha 
achievement of such “section aeroisoclinicism" is very dif- 
cult and can be realized only by separate variation of the 
bending and torsion stiffnesses; even so, the aeroisoclinic 
condition obtains for only one type of aerodynamic loading 
condition at one Mach number. However, an overall type 
of aeroisoclinicism in which bending and torsion action tend 
to cancel for the wing as a whole is relatively easy to achieve. 
This overall type has, for practical purposes, the same 
advantages аз section aeroisoclinicism, in that the аего- 
elasti¢ phenomena considered in this report tend to be 
negligibly small for such a wing. 

As may be seen from figure 2, at a small positive value of 
the parameter k the values of the parameter 4% for diver- 
gence as well as those for given shifts in the spanwise center 
of pressure tend to infinity. This particular value of Е 
represents aeroisoclinic wings in the overall sense; from 
equations (18) and (19) 15 тау be seen to be the reciprocal 
of the value of Ка given ір table 2. Hence, from the definition 
of k (eq. (9)), - 

8: (GJ), 1 


€1C; (ЕТ), inn M 


with the implication that the distributions of the stiffness 
are of either Ше с or the constant-stress type and that К. 
pertains to either of these types and to the appropriate taper 


(33) 


. ratio. Equation (33) indicates that, for а given plan form 


with assigned values of 8;, с, and A, the disposable param- 
eters for the achievement of aeroisoclinicism are the elastic- 
axis location e, which enters into the parameter e;, and the 


(GJ), 


root-stifiness ratio (ЕТ) ; the aerodynamic center is not under 
Т 


the control of the designer to any appreciable extent. 

А decrease in the torsional stiffness can sometimes Бе 
effected without decrease in the bending stiffness or impair- 
ment of the strength characteristics of the wing, and overall 
aeroisoclinicism may be achieved in this manner for swept- 
back wings. Or, if aerdisoclinic conditions are considered 
at the outset, a wing can be designed with the elastic-axis 
location relatively far back (in the case of & sweptback wing) 
or forward (in front of the aerodynamic center, in the case 
of a sweptforward wing) in order to achieve aeroisoclinicism. 
However, the fact that only certain types of aeroelastic 
phenomena are considered in this report must be kept in 
mind. Locating the elastic axis far back or decreasing thc 
torsional stiffness, for instance, may lead to flutter difficul- 
ties, the solution of which may require excessive mass 
balancing of the wing as a whole. 


RELATION OF THE CHARTS TO DESIGN PROCEDURE 


The first step in the design of a wing structure, once the 
wing geometry and the, overall airplane characteristics 
have been decided upon, usually consists of a rough appor- 
tioning of structural material along the span in a manner 
intended to satisfy strength requirements approximately. 
At a later stage in the design procedure the structure it 
checked for aeroelastic effects and modified, if necessary. 
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The modifications are then checked again, and so on, until 
both stiffness and strength requirements are met with what 
is believed to be а near-optimum structure from weight 
considerations. The charts of this report may be used to 
facilitate the procedure at several stages. 

At the very outset, the preliminary-survey charts can be 
used to establish some overall aeroelastic characteristics of 
the wing structure that would be obtained by designing the 
wing for strength alone. If these characteristics are satis- 
factory, the design of the wing structure can proceed on the 
basis of strength requirements alone. The final design can 
then be checked for the aeroelastic effects considered in this 
report by means of the charts contained herein, and for other 
aeroelastic effects, such as flutter and loss of lateral control, 
by equally approximate methods. However, if the pre- 
liminary survey indicates that a wing designed on the basis 
of strength alone would be unsatisfactory from consideration 
of аегойавис ву, sufficient additional stiffness may be incor- 
porated in the preliminary design stage, provided the taper 
ratio does not differ greatly from 0.5 and the wing thickness 


ratioisconstantalong thespan. Forinstance, the preliminary- | 


survey charts may indicate a shift in the spanwise 
center of pressure which gives rise to a shift of 4 percent їп 
the aerodynamic center; whereas the desired maximum shift 
is 2 percent, so that the spanwise shift must be reduced 10 
50 percent of that indicated on the preliminary-survey chart. 
The shifts in the spanwise center of pressure for a wing with 
increased stiffness at the tip (the “excess strength” case, for 
which о--2,0) and for a wing with decreased stiffness at the 
tip (the wing with w=0.5, for which the results of the case of 


Ыр = 0.5 may be used) сап then be obtained from figure 


9 (b) and equation (32), in conjunction with the value of the 
dynamic pressure at divergence estimated from equation (18) 
ог (19). The fact that the wings with о--2.0 and w=0.5 
have different dynamic pressures at divergence than does 
the constant-stress wing must be kept in mind. 

From the shifts of the spanwise center of pressure cal- 
culated ш this manner the value of о for the desired spanwise 
ghift can be obtained by interpolation and, hence, the approx- 
imate magnification factors to be applied to the stiffness 
distribution for constant stress can be obtaimed from the 
chart at the top of figure 14. Estimates for the other 
aeroelastic characteristics considered in this report can then 
be obtained for the wing with this modified stiffness distribu- 
tion by interpolating between the results given for these 
characteristics for wings with w— 0.5, 1.0, and 2.0; that 18, 
for the cases referred to, respectively, as 


(М с); m 

uo e 

(Ме. 
ie m 


апа 


Come (excess strength) 
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in table 2 and figures 5 (b), 7 (b), 9 (b), and 11 (b). Once 
the structure of such & wing has been designed, the various 
aeroelastic effects considered herein should be checked by & 
more accurate method, such as that of reference 1, and the 
loss of lateral control and the flutter characteristics should 
be calculated. 


CONCLUDING REMARKS 


Charts have been presented for the estimation of аего- 
elastic effects on the spanwise lift distribution, lift-curve 
slope, aerodynamic center, and damping in roll of swept 
and unswept wings at subsonic and supersonic speeds. 
Two types of stiffness distributions have been considered: 
one which consists of a variation of the stiffness with the 
fourth power of the chord and is appropriate for solid wings, 
and one which is based on an idealized constant-stress struc- 
ture and is believed to be more nearly representative of 
actual structures. 

The limitations of these charts are that they do not apply 
to wings with very low aspect ratio or very large angles of 
sweep nor to wings with large sources of concentrated аего- 
dynamic forces. The charts are likely to be less reliable 
for wings with zero taper ratio than for wings with other 
taper ratios and less reliable when the component of the 
Mach number perpendicular to the leading edge is transonic 
than when this component is either subsonic or supersonic. 
Wings with large discontinuities in the spanwise distribution 
of the bending or torsional stiffnesses cannot be analyzed 
directly by use of the charts, but a means of making approxi- 
mate calculations for such wings has been presented. No 
charts have been presented for inertia effects but a method 
of estimating these effects has been outlined. 

In addition to facilitating the calculation of various static 
aeroelastic phenomena, the charts serve to simplify design 
procedure in many instances, because they can be used 
at the preliminary design stage to estimate the amount 
of additional material required to stiffen a wing which is 
strong enough and because they indicate that the best way 
of distributing this additional material is to locate most 
of it near the wing tip. 

Also, the charts facilitate the achievement of aeroiso- 
clinic conditions, inasmuch ая they serve to define a simple 
relation between the elastic-axis location and the wing 
stiffness ratio which is required to obtain this condition 
for a given plan form. Finally, the charts indicate that a 
wing which is strong enough is most likely to be affected by 
aeroelastic phenomena if it is to operate at high dynamic 
pressures, if it is thin, if it has a large angle of sweep, if it is 
designed for a low wing loading, if it has an elastic-axis 
location relatively far back on the chord, and if it is to 
operate at transonic or high supersonic Mach numbers. 
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METHODS OF CALCULATIONS ON WHICH THE CHARTS ARE BASED 


THE AEROELASTIC EQUATIONS 


The methods of calculating aeroelastic phenomena used in 
preparing the charts of this report are based on the following 
assumptions: - 

(1) Aerodynamic induction is taken into account by apply- 
ing an overall correction to strip theory and, when matrix 
integrations &re used, by rounding off the resulting load 
distribution at the tip. 

(2) Aerodynamic and elastic forces are based upon the 
assumption of small deflections. 


(3) The wing is clamped at the root perpendicular to a 


straight elastic axis (see fig. 1), and all deformations are con- 
sidered to be given by the elementary theories of bending 
and torsion about an elastic axis. 

In keeping with assumptions (1) and (2), the force per unit 
width on & wing section perpendicular to the elastic axis is 


a7 (Сг, æt Cr, 4) (A1) 


а 
where a, and a, are, respectively, the angle of attack due to 
structural deformations and the rigid-wing angle of attack, 
in planes parallel to the plane of symmetry. (The geometri- 
cal angle of attack is considered to be constant along the span 
in equation (À1); in the case of linear twist the coefficient 
Or. is used instead of Съ.) The torque of this force about 


the elastic axis is Lec for uncambered sections. 
The integral equations for the accumulated torque and 
the bending moment are 


r-[" lec ds 


м-| "(tds ds 


(A2) 


(Аз) 


and, insofar ав assumption (3) holds, Ше angles of structural 
twist and bending referred to axes parallel or perpendicular 
to the elastic axis are 


=f P (A4) 
41 
r-f gr Мав (A5) 


The angle of attack due to structural deformations is related 
to e and Г by the equation 


аұс-ф cos A—T sin А (A6) 
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Combining equations (А1) to (А6) gives two simultaneous 
differential equations: 


— qec? 


LAL [Cz ast Со бе сов A—T sin AJ] (A7) 


s (07 5 - 


та (2192 25 щил "n 1 = [Cz % + Cr, КС cos А—Гзш А)| (А8) 


These equations are subject to the following boundary 
conditions: 
Zero twist and bending at the root, 


Ф(0)-+0 (Ада) 
T(0)—0 (А9Ъ) 
Zero torque, moment, and shear at the tip, 
de 
(47 unt (A102) 
( EI dr —0 (A10b) 
ds ses r 
dir 
(z1 9) _,=0 (А10с) 


In the following sections, equations (А7) апа (А8) are 
solved explicitly for ап untapered wing with constant stiff- 
ness along its span and by matrix integration for а wing with 
any arbitrary stiffness and chord variation. 

SOLUTIONS FOR UNIFORM WINGS 


Арыс geometric angle of attack.—If the torsional 
stiffness, the bending stiffness, and the chord of the wing 
have constant values (GJ),, (ЕГ),, and с,, respectively, along 
the wing span, equations (А7) and (А8) become 


p” cog A— — q* E с + (e cos А-Г sin М (А1) 


Y din дез E (оова Ра М (А12) 


where the differentiation denoted by the primes is with 
respect to £z 1—5 and the dimensionless parameters 0% and 
4 are defined by 
| Съ 4616,28,8 cos А 
D— Tia GIy 
Си 96:88 sin А 


а (ED, Н (414) 


(A13) 
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Differentiating equation (А11) onco with respect to Е and 
combining it with equation (412) yields the single differential 
equation 

, 


/ 
o, 4- g* o, — do, — q* 2+9 = (А15) 


(Тһе factor 2 is used with a, for the sake of consistency, 


despite the fact that a geometrical angle of attack which is 
constant over the span does not have a spanwise derivative.) 
Equation (А15) is subject to the following boundary 
conditions: 

From equations (Ада) and (A9b) 


os (1)--0 (A16) 
From equations (A102) and (A10b) 
os! (0) —0 (А17) 
From equations (À10c) апа (А11) 
a!" ()=—g* | 4-0, (0) | (A18) 


Where functional notation is used, so that, for instance, o, (1) 
means the value of o, at £—1. 


The solution of equation (А15) can be effected very readily | 


by means of Laplace transforms. The complete solution of 
this equation is 


H (1) 


о =F) А(9-Н(9 | ш 


where Ше integral H(£) is defined ав 
| 
HOE f | AAE- -TEE (420) 
The functions fs, fı, and fs are defined by 


Ja (£) = Cie ев (0. сов vet sin yt) (A212) 


Ja) = Оце № eft (о, сов ү+@ віп yt)  (A21b) 


fs (B = Сет? рев С соз veto sin vE) (А21с) 


where —28 and В-Ту are Ше roots of 7*--q*r—g—0 and 


__ 482 
ба, r3 
O= 1 -(, 
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The condition {от divergence is that о, be finite when a, is 
zero along the entire span. Ав can be seen from equations: 
(А19) and (À20), divergence can occur only when 


Ж(1)--0 


Thus, for а particular value of the parameter к=, the value 


(A22) 


044% (or 4) at divergence is the one which satisfies equation 
(A22). 

Constant geometric angle of attack.—For the additional- 
type angle of attack, о;(Е)-- Constant: 


Н @=[1— «6 + ‚ (A23) 


and 


ков (Е) ох fE) 
Og 21) 


For lift distributions based on assumption (1) given at the 
beginning of this appendix, the lift per unit width of span 
m&y then be written as 


(A24) 


1 fW 
kh f) 


The wing lift coefficient, the wing-root bending-moment 
coefficient, and the wing-root twisting-moment coefficient are 
given in general by 


Le 


бы== FS 


= Cr, в ССС 


0 


---- — + 
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4M, 
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S es J, ], E iro eot ae 


== Оу, (А2 6) 


81€; 


= 0r, 35а ), (2) тка ола 


Then, for the uniform wing, 
0 Oy T 


АС) 
~ FD) 


Ge? ы eae ae 
B 0/0 Og 


— o а) 
2 S) 


On _ 
Cr, 


(A27) 


(A28) 


(A29) 
and | 

Or. 
Ch. 


linearly varying geometric angle of attack.—For Ше. 


linear-twist-type angle of attack, «,(f)=(1—é)a,,, the factor 


xis 1, and 
H(i)=a,()— [309 —714(8] о, (A30) 
so that 
eas al +оо 
я о, 
(1 
= Ә-Ә (A81) 


'The ratios of the wing lift, wing-root twisting moment, and 
wing-root bending-moment coefficients to their rigid-wing 
values are then, on the basis of assumption (1), 


Cr. Кер 5 чии 
ne 2 |, wD Га) ДЕ 


x ES ні (А32) 
0-48 CAE) 
E= hh ae dt 
EE =з} т“ (1) + 9701) — ці (А33) 
апа 
Ст Зара Cr, 
С. б 


ав in the preceding section. 
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SOLUTIONS FOR NONUNIFORM WINGS 
Equation (À1) may be written in the matrix notation of 
reference 1 as 


COTA Ora, 


144 (434) 


afes 


and equations (А2) and (А3) as ши 
(T) еве; 


= hp Cr ео (E) | fiait ta} (439 
(M ms ATA) 


сев zelo] аз 
where the matrix [//] performs an integration of the running 
torque [ес from the tip inboard, and the matrix |7171 per- 
forms & double integration of the running load from the tip 
inboard. ‘These matrices are derived and given in reference 
1. They are based upon Simpson’s rule with а modification 
at the tip, where the load distribution is assumed to go to 
zero With an infinite slope at the tip. 

Equations (A4) and (A5), written in matrix notation, are 


ur GF | cn) 


a 


д” 


(T) ЕЕ; (A38) 


where the matrix [1]”’ serves to integrate the accumulated 
torque or bending moment outboard from the wing root. 
This integrating matrix is based upon Simpson's rule without 
the tip modification and is given in reference 1. . 

The substitution of equations (A35), (А36), (A37), and 
(438) in the matrix equivalent of equation (А6) yields 


(а) = 14) Сои) а ^ (830) 


where the aeroelastic matrix [4] is defined by 


ЕЛІНЕН 


The parameters 0% and 4 are denned by equations (A13) and 
(A14), respectively, and 


к= 
When a, is zero along the entire span, equation (A39) 


becomes | 
{в,}=@*[А]{о,} 


Consequently, for а particular value of b, the yalue of 4% 
&t divergence can be found by the iteration of the aeroelastic 
matrix [4]. 


(A41) 
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Equation (А39) may be rearranged as follows: 
[[1]— g*[4]] («es - o4] (ак! 


The set of linear simultaneous equations represented by 
equation (А42) may then be solved for the total angle of 
attack ко, |- o, in terms of the values of a, along the span. 
The integrations in equations (А25), (А26), and (А27) 
тау be performed with the first rows of the [77] and [777] 


matrices. Thus 
бы h| E irata) 


(A42) 


ta h| È Jta "= 
л E (177 is | коса o } dus 
го КЛЕТ 
ала 
Pale cone 


The aeroelastic characteristics of uniform wings were 
calculated by both the direct method of the preceding section 
and the matrix method given in this section. ‘The values of 
the divergence parameter 4%р, calculated by the direct 
method, were found to be about 5 percent greater than the 
corresponding values calculated by the matrix method. 
This discrepancy can be shown to be almost entirely due to 
the rounding off of the loading of the wing tip in the matrix 
method. The differences between corresponding values of 
©; Съ Ов Ст ae 
a Cin’ д. and Cr, are negligible. 


0 


COMBINATION OF RESULTS 


The forms of the approximate formulas used in combining 
the results of the many computations indicated in the analysis 
were obtained by considering a highly idealized semirigid 
wing; that is, a wing which is rigid along its entire span but 
can bend and twist at the wing root subject to the restraint of 
а bending and а torsion spring. 

If it is assumed that the two spring constants correspond 
to (64); апа em, the value of g* at divergence is given by 

t { 
the simple formula 
K. 
* __ I 
даа (A40) 
where the factors A, and Ка depend on the taper ratio and 
the spanwise variation of the stiffness. As shown in reference 
2, this formula serves as a good approximation to the calcu- 
lated values of q*». 


321005—55———00 
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For the semirigid wing, the ratio of a, to a, is found to be 


proportional to Янг. In order to adapt this expression 
ix 


Фр 
to Ше flexible wings considered in the present analysis, the 
following approximate expression was found to provide 
satisfactory correlation: 
i. 
Ses HP. (ГРАД) (A47) 


а, ку 


9р 


where f and Af are functions of the spanwise coordinate 8% 
and the wing-chord and stiffness variations; F is a function 
of the parameter k and the wing-chord and stiffness varia- 
tions. The functions f, Af, and Е also depend on the type 
of spanwise variation of the geometrical angle of attack, the 
subscripts 1 and 2 being used to distinguish between the two 
types of interest. Тһе accuracy of equation (А47) is 
illustrated in figures 16 and 17. 

If equation (А25) is used for the wing lift coefficient (with 
Е replaced by s*) and equation (А47) for the angle-of-attack 
distribution, the wing lift coefficient may be expressed as 


q 
С, шеру gp C. +O, (А48) 
Яр 
or 
q 
Cr, aa (1--») 
бо, mi (A49) 
dp 
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Figure 16.—Comparison of angle-of-attack ratios calculated by the 
matrix method of appendix А with those calculated from equation 
(20) for constant geometric angles of attack at various dynamio- 
pressure ratios. Х=0.5; Е--8; stiffnesses proportional to с“. 
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Етаувв 17.—Comparison of angle-of-attack ratios calculated by the 
matrix method of appendix A with those calculated from equation 
(20) for constant geometric angles of attack and for various values 


of the parameter Б. Х--0.5; i1 stiffnesses proportional to c. 


where the parameter 


i (f+F Af) ds* 
мы ее (А50) 
| ге Og ds* 
о Cr 477 
Cr | : 
во that » and G = are functions of Ё and of Ше wing-chord and 
L 


wo | 
stiffness variations and depend on the type of geometrical 
angle-of-attack distribution as well. 


As indicated by equation (А47), within the approximation | 


inherent in that equation, the shape of the spanwise distribu- 
tion of а, does not vary with dynamic pressure. Therefore, 
to a good approximation, the lateral center of pressure of the 
lift due to o, (as well as that due to œz) does not change its 
position along the elastic axis when the dynamic pressure 
changes. The following approximate formula for the wing- 
root bending-moment coefficient тау therefore be deduced 
from equation (А48): 


4 | 
2 Съ HE Org, 


А51 
и (А51) 
до 


Cp=8,*v 


where. 8,* and 8,* are the dimensionless moment arms about 
the effective wing root of the lifts due to a, and a, and are 
defined by : 

M 


8, = 


1 
9 Ly 81 


and 
M, 
8, = 1 Е 
2 Г. в, 
ТЪеп, 
1—-—(1—u») 
С, 4 А 
—— = 2 
Co, A (А52) 
dp 
where и is defined by 
8," 


во that Са. is а function of the parameter Ё, of the taper ratio, 
B, 


and of the stiffness distributions; it also depends on the type 
of geometrical angle-of-attack distribution. 


Ап approximate . formula for the wing-root twisting- 
moment coefficient may be deduced from equation (А48) as 
follows: 


а Or (А54) 


where e, and е , are Ше effective dimensionless moment 
arms &bout the elastic axis of the lifts due to o, and o, and 
are defined by 


and 


Then 
ит Ч (A65) 


where 
тез 14 (А56) 


so that oe is а function of Ё, the taper ratio, and the stiffness 
Т, 


0 
distributions and also depends on Ше type of geometrical 
angle-of-attack distribution. 
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The values of р, и, and т are given for the two types of 


geometrical angle-of-attack distributions in figures 9 and 11. 


Figure 18 shows the approximate formulas (А49), (A52), 
and (А55) to be in good agreement with more accurately 
computed values. | | 

The foregoing approximate formulas for the structural 
angle of attack and for the lift, bending-moment, and 
twisting-moment coefficients are not applicable to wings with 
zero taper ratio. Ап attempt was made to combine and to 
systematize the results calculated. for such wings in the 
manner employed for wings with other taper ratios, but the 
approximate formulas obtained in this way were found to 
yield unreliable results. Consequently they аге not pre- 
sented in this report; instead, the results calculated for Ше 
wings with zero taper ratio are presented directly in figures 
6, 8, 10, and 12. 
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Figure 18.—Comparison of lift- and moment-coefficient ratios for 
constant geometric angles of attack calculated by the matrix method 
of appendix A with those calculated from equations (26), (27), and 
(28). А--0.5; k=8; stiffnesses proportional to сі. 


APPENDIX B 


STIFFNESS DISTRIBUTION OF CONSTANT-STRESS WINGS 


OUTLINE OF CONSTANT-STRESS CONCEPT 


In order to calculate aeroelastic effects, the bending and 
torsional stiffnesses of the wing structure ЕТ and GJ have to 
be known. These stiffnesses enter the calculations in two 
ways. The root stifinesses, as indices of the overall bending 
and torsional stiffnesses, constitute primary parameters 
which are required in the use of the charts of this report but 
were not required in the preparation of the charts. On the 
other hand, the stiffness distributions—that is, the ratios of 
the local stiffnesses along the span to the root stiffnesses—are 
secondary parameters which are not required in detail in the 
use of the charts but did have to be assumed in order to pre- 
pare them. 

In calculations preliminary to the actual design of the 
structure, the bending and torsional stiffnesses of the struc- 
ture are not known; they must be estimated on the basis of 
either past experience or considerations of an idealized struc- 
ture. For the purpose of estimating stiffness distributions, 
past experience with similar structures is likely to be а useful 
guide in any specific case but does not lend itself to generali- 
z&tion апа hence to the preparation of generally applicable 
charts. The stiffness distributions (other than those which 
vary as the fourth power of the chord) used to prepare the 
charts of this report have been obtained from considerations 
of an idealized structure, as outlined in this appendix. 

Basically, the method of this appendix consists in an effort 
to relate the stiffness of a wing to its strength and to estimate 
that strength on the basis of certain assumptions. The 
fundamental assumptions are that the bending and torsional 
stresses аге constant along the span and that the applied 
loading is proportional to thelocalchord. The other assump- 
tions concern the bending and torsional stresses caused by 
this load and their relation to their allowable values. Іп 
estimating these stresses the structure is assumed to be 
essentially of the thin-skin, stringer-reinforced shell type. 
Certain effectiveness factors are used—for instance, the ratio 
of the allowable torsional stress to the allowable bending 
stress, or the ratio of the cross-sectional area of the effective 
torsion cell to the product of the chord and the wing thick- 
ness. The root stiffnesses estimated by the method of this 
appendix depend directly on the values of these ratios. The 
stiffness distributions, on the other hand, are largely inde- 
pendent of these ratios but imply the assumption that the 
ratios are approximately constant along the span. Conse- 
quently, the constant-stress concept used in this appendix is 
more likely to furnish useful results for stiffness distributions 
than for the root stiffnesses, and, because of the type of 
structure assumed, the concept is not applicable to very thin 
wings. 
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ASSUMED APPLIED LOADS 


If the applied normal load is distributed in a manner 
proportional to the chord, that is, if 
l= Ke | - (B1) 
the bending moment at any point on the span can be obtained 
by integrating the chord distribution asfollows: 


M=K аур c ds* ds* 


where 8% is the dimensionless distance along the reference 
axis measured from the effective root. Similarly, the total 
normal load on one wing is given by 


g X12 " 
-кав | сав 


If the wing is linearly tapered, so that 
с=с, [1—(1—)в*] 


(B2) 
where the taper ratio Х is defined by 


y set 
Ст 
then the ratio of the bending moment at any point of the 
span to the product of the total normal load and the wing 
semispan less one-half of the fuselage width can be expressed 
as follows: 


рве) (B3) 


cos À 


where the function f, of 8% and А is defined by 
14-24 1-2 st) (1—в*): 


8 IA 14А 


(B4) 
and shown in figure 19. 

Similarly, if the moment arm of the normal Joad applied to 
the wing аф any station is also proportional to the chord, the 
constant of proportionality being e;, the distributed torque 
at any station is then e,cl, and the accumulated torque is 


Тт-ак 22. ( ади (В5) 
which may, in turn, Бе expressed as 
Т, 
ЕР" Л (8*, X) (B6) 
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where the function ў; of 8% and А is defined by 


APIA (2-5 0—» a N ры 
Шегі Кгееу арии * теу 19-9) (ВО 


апа the average chord с is defined by 


= (В8) 


The function fy is shown in figure 20. 


The total normal load on one wing P can be estimated 


from the design gross weight and the design load factor of the - 


airplane in the following manner: 


If the fraction of the wing lift to the total lift carried by the 
airplane (including that of the fuselage and tail) is designated 
by 71, ВО that 
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Етсовв 19.—The bending-moment function fs. 
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and the fraction of the wing weight (including the amount of 
fuel, external stores, and so on used in the critical design 
condition) to the total design gross weight is designated by 
12, 


tb 


Mi 


72— W 


then equation (B9) may be written as 


Pas "i (B10) 


where 
13—71: 18 


With the value of P given by equation (B10), equations (B3) 
and (B6) serve to express the local bending and torsional 
moments in terms of known design parameters. 


EFFECTIVE SKIN THICKNESS REQUIRED TO RESIST APPLIED LOADS 


The wing structure has to resist both the applied bending 
moments and the applied torques; in other words, the load- 
carrying members must resist combined axial and shear 
stresses. A relation commonly used in the design of wing 
structures loaded by compressive and shear stresses due to 
bending and torsion moments is 


TOR 


where fo is the applied bending stress, f: the applied shear 
stress, Ев the allowable (compressive) bending stress, and F7 
the allowable shear stress. However, & similar relation, 


f GE 
Re (В11) 


is somewhat conservative and much more convenient for Ше 
present purpose апа, consequently, is used as the basis of the 
following development. If the margin of safety is not zero, 
equation (B11) can be rewritten as 


Pitan (B12) 


where 7, is an effectiveness factor which can be expressed in 
terms of the margin of safety (М. 5.) аз 


1 


ТЕМ. 8 ы, 
ТЪе applied bending stress is 
ре (B14) 


where z is the maximum ordinate on the compression side 
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measured from and normal to the chordwise principal axis. 
Similarly, the applied shear stress is 
Т 

f: 228 (B15) 
where A is the cross-sectional area of the (assumed) single 
torsion cell and # is the skin thickness on the compression 
side. Substitution of equations (B14) and (B15) into 
equation (B12) yields 


. Mz 
IF, (1+ 


T Е I 


M Fe 22u) " г 


In order to relate the bending and torsion stiffness of the 
wing to the skin thickness і or to an equivalent thickness & 
which includes the material in the stringers and spar flanges, 
the bending stress 15 assumed to be carried by а, box covered 
with sheet of an effective thickness &, the webs are assumed 
to carry no bending stress, and the torsion stress is assumed 
to be resisted by an equivalent single cell, the two webs of 
which contain all the material of the actual webs. The tor- 
sion and bending moments of inertia may then be written as 


2-4 At 
р 7 9 
2 
=92 1817 ору, 17 
7879 T (В 
апа 
hN? h 
I= тис! Е (5) nu? + T1213 14? ($) | 
БАС | 
— 7156 (5) t (B18) 


where the effectiveness factors 7; to 715 are defined in table 
1. Іп the factor 1, the effective perimeter 2 of the torsion 
cell is the sum of the lengths of skin around the perimeter, 
each weighted by the ratio of the thickness of the critically 
stressed element to the thickness of the given length of skin. 

When the value of I given by equation (B18) is substi- 
tuted into equation (B16), equation (B16) may be written as 


T Fa xO (B19) 


t= 
UM Fr 242 


Mz | 
БИШНЕ” ия 
| Ғ BN4715C (5) 


By making use of equations (B3), (B5), and (В10), as well 
as the effectiveness factors 716 to тъ defined in table 1, 
equation (B19) can be written as 


nans nW 612 fe (1 тава 1 
өүө в c? соз А с д, FA een 


The factor f; is defined in terms of the factors f, and f, given 
by equations (B4) and (B7) as 


(B21) 


The function ўз is shown in figure 21. 
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Еісоке 21.—The moment-ratio function fs. 


BENDING AND TORSIONAL STIFFNESSES 


Substitution of the value of £, given by equation (B20) 
into equation (B18) yields an expression for the bending 
moment of inertia Г or for the bending stiffness ZI at any 
point along the span. The value of this stiffness at any 
point on the span may be divided by the stiffness at the 
wing root (EI), This ratio can then be expressed ая 


T 
вт вала A. 
(ЕГ), DEDE n 


__ Ме 7 1961 
=: (д, HE) 928 
where 
f 7 1961 
8 

pace 2-4 (В23) 

Cr fr, 71961 

Limo 


The function f»isplotted in figure 22. The value (ZI), may 
be obtained from equations (B18) and (В20) 88 


(ED,= та É АУ ы О в ие! (В24) 


where 


i (1A) fo, 7196 
га (нац) 


(В25) 
The function Е, is shown in figure 3 as а function of À with 
719€1 
Аһ 

Similarly the torsional stiffness GJ may be obtained by 
substituting the value of ¢, given by equation (B20) into 
equation (B17). However, from equations (B17) and (B18) 
the ratio of the torsional stiffness to the bending stiffness 
may be obtained in the form 


GJ ай nent 
EI ТЕ asnonis 


ав the parameter. 


(В26) 
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This equation shows that the ratio GJ/EI is constant along 
the span within the framework of the constant-stress con- 
cept. Equation (B26) may, therefore, be interpreted as an 
expression for the value GJ'/EI at the wing root, that is, for 
the value (4J), (EI), The torsional stiffness at any other 
point on the span can then be obtained from equation (В22), 
since 
GJ | EI 
(GJ). (ЕТ), 


because GJ/EI is constant over the span. 

The stiffness ratios GJ/(GJ), and EI/(EI), can be ob- 
tained directly from figure 22 when the thickness ratio Л/с 
of the wing is constant along the span; if the thickness ratio 
is not constant, the factor f, obtained from figure 22 must 


be multiplied by the ratio n J at any station to obtain the 


stiffness ratio at that station. As may be seen from figure 
22, the function fə does not vary much with the parameter 


285, this parameter represents the additional amount of 


(827) 


skin thickness required to carry the torque (see eqs. (B19) 
and (B20)), and this additional thickness is small for most 
conventional wing structures. Consequently, an average 


value of Ln £1971.:0.03 was used to obtain the stiffness distri- 


butions used in the aeroelastic calculations on which the 
charts of this report are based. 
Equation (B22) shows that, once & value has been assumed 


for the term таа, the stiffness ratios EJ/(EI), апа GJJ(GJ), 


gre independent of the effectiveness factors used in this 
analysis. "Therefore, specific values of these parameters need 
not be known in order to estimate the stiffness distributions, 
but one of the assumptions on which equation (B22) 1s based 
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is that whatever values the effectiveness factors have are 
nearly constant along the span. In order to estimate the 
value of (JZI),, however, these factors must be known, since 
they enter directly into equation (B24). Estimates of (ЕТ), 
and (6.7), obtained in this manner are, therefore, subject to 
all the limitations imposed by the approximations.of the 
constant-stress concept. Hence, some judgment must be 
exercised in using these estimates, end, if possible, they should 
be modified in the light of experience. 
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STRUCTURAL WEIGHT ASSOCIATED WITH THE STIFFNESS DISTRIBUTION 


The increase in structural weight associated with a given 
increase in stiffness can be estimated on the basis of assump- 
tions similar фо those made in relating the stiffness and the 
strength. For Ше purpose of this analysis the various com- ` 
ponents of the wing structure, exclusive of the carry-through 
structure within the fuselage, are classified in two groups: 
one which contains the elements that t&ke the bending and 
torsional loads due to the assumed loading and one which 
contains all other components. Іп the first group are 

(1) The amount of top and bottom skin that is used in the 
estimation of the thicknesses required to withstand the 
bending and torsional loads, including stringers and spar 
flanges included in the equivalent skin 

(2) Webs, including any web stiffeners 
In the second group are 

(1) Skin, stiffeners, false spars, and so on, which are not 
considered in the estimation of the equivalent thicknesses 

(2) Ribs, bulkheads, and posts designed to raise the buck- 
ling strength of the cover sheets 

(3) Control surfaces and their supports, attachments, 
and actuating mechanisms 

(4) The supports of internal по 

This analysis is concerned only with Ше first group and, 
more specifically, with the relative increase in the weight of 
this group occasioned by an increase in stiffness of the main 
structure. Means of estimating the actual magnitude of the 
weights involved and of estimating the weights of some of the 
items in the second group as well are given in references 9 
and 10. 

The weight per unit ден of the structural elements of Ше 


` first group can be written as 


(В28) 


where y, is the density of the material of the primary struc- 
ture (or ап equivalent density in the case of sandwich con- 
struction), па is the ratio of an equivalent perimeter 2 to 
the actual perimeter of the cell, and 2 is the sum of all the 
lengths which constitute the perimeter, each multiplied by 
the ratio of its equivalent thickness to the. equivalent 
thickness 1, of the upper cover sheet. | 

In view of the assumption made concerning the combina- 
tion of bending and torsional stresses, the thickness i, re- 
quired in equation (B28) can be obtained from equation 


= 290872146 


` (B18) as 
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Consequently, the total weight (for both wings) of the struc- 
tural elements of the first group can be estimated from the 


relation 
W, 2 = 1 І IH. fe E 
"WI Js (ЧЫЗ (529: 
* cos А . 


Equation (B30) serves to estimate relative changes іп Ше 
Weight of the first group of structural elements. For instance, 
with a given distribution of J and Л, that weight is directly 
proportional to J, and inversely proportional to 22. Simi- 
larly, given two different distributions of J and À with the 
same values at the root, the ratio of the weights is equal to 
the ratio of the two values obtained by using the respective 
distributions of J and h in the integral of equation (B30). 

Although the actual value of W, is not relevant to this dis- 
cussion, it may be estimated by substituting the previously 
calculated stiffness distributions into equation (B30), and 
the result is given here as a matter of general interest: 


"атыла Ys nW 57/2 
uc 74715718 Ев dye Ах А соз ДО” (B31) 
Cr 
where 
ГИ. Gok 
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According to equation (B31), the structural weight is di- 
rectly proportional to the design gross weight, load factor, 


REPORT 1140—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


swept-span aspect ratio, span, and density of the material 
of the primary structure and inversely proportional to the 
allowable stress and the wing thickness ratio. The depend- 
ence of the weight on the taper ratio (all other parameters, 
notably the aspect ratio and span, are the same) is illustrated 
in figure 14 by a plot of the function Е against taper ratio for 
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several values of Ше parameter я апа for several ratios of 


the wing thickness ratios at the tip and at the root. 
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